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We describe a method for determining y using B decays followed by a multibody D decay. In the talk we focus on K$ K~ K + 

final state, but other modes such as D — > K$K~K + and D — > K$ iz~ft + iz° can also be used. The main advantages of the method are that 
it uses only Cabibbo allowed D decays, and that large strong phases are expected due to the presence of resonances. Since no knowledge 
about the resonance structure is needed, yean be extracted without any hadronic uncertainty. 



1 Basic idea 

The theoretically cleanest way of determining the angle 

Y=M-VudV^/V cd V* b ), (1) 

is to utilize the interference between the b — > ciis and 
b — > ucs decay amplitudes [1-7], see Fig. 1. The salient 
feature of these transitions is that they involve only distinct 
quark flavors and therefore do not receive any penguin con- 
tributions. In the original idea by Gronau and Wyler (GW) 
[1] the — > DqpK^ decay modes are used, where Dqp 
represents a D meson which decays into a CP eigenstate. 
The dependence on y arises from the interference between 
the B ± -> D ^ and B ± -> 25° A"* decay amplitudes. The 
main advantage of the GW method is that, in principle, the 
hadronic parameters can be cleanly extracted from data, by 
measuring the B^ — > D°A ± and B^ — > D°A ± decay rates. 




Figure 1. The dependence on / arises from the interference between the 
B ± -> D ^ and B ± -> TjV* decay amplitudes. 

In practice, however, measuring y in this way is not an easy 
task. Due to the values of the CKM coefficients and color 
suppression, the ratio between the two interfering ampli- 
tudes, rs [see Eq. (4)], is expected to be small, of order 
10% — 20%. This reduces the sensitivity to y, which is 
roughly proportional to the magnitude of the smaller am- 
plitude. It also leads to experimental difficulties in measur- 
ing the color suppressed B — > D°K~ mode (and its charge 
conjugate) preventing a straightforward application of the 

"Talk given by J. Zupan, based on [8]. 



GW method [2]. There exist a number of extensions of the 
original GW proposal which avoid the problem by not rely- 
ing on the measurement of B -^D°K~ amplitude [2,3,6]. 
Instead several different decay modes of D mesons such 
as quasi two-body D decays with one particle a resonance 
(e.g. D° — > K* + n~ [2,3]) are used. Since these are really 
three body decays (for instance in the example mentioned, 
K* + decays strongly to K°n + or K + n°) one can pose the 
following questions : 

• Can one use the complete phase space of such three- 
body D decays for y extraction? 

• Is it possible to avoid fits to Breit-Wigner forms in 
doing the Dalitz plot analysis? 

As we show in the following, the answer to both of these 
questions is positive. The first question was raised already 
in [3], however, most of the results and applications we 
present are new. For the sake of concreteness, we con- 
centrate on the D — ► K$n~n + decay mode, while an ex- 
tension to a larger set of the decay modes can be found 
in [8]. The advantage of using the chosen decay chain is 
threefold. First, one expects large strong phases due to the 
presence of resonances. Second, only Cabibbo allowed D 
decay modes are needed. Third, the final state involves 
only charged particles, which have a higher reconstruction 
efficiency and lower background than neutrals. The price 
one has to pay is that a Dalitz plot analysis of the data is 
needed. We describe how to do the Dalitz plot analysis 
in a model-independent way, and explore the advantages 
gained by introducing verifiable model-dependence. The 
final balance between the advantages and disadvantages 
depends on yet-to-be-determined hadronic parameters and 
experimental considerations. Finally, we mention that an 
equivalent formalism to the one we present below has been 
independently developed by Atwood and Soni in [9]. 
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2 Model independent determination of y 



Let us focus on the following cascade decay 2 

B ^DK -> (K s n~n + ) D K- , 

and define the amplitudes 

A(B- ->D°K-) = A B , 

A(B- -» dV") ee A B r B e i{SB -V\ 



(2) 

(3) 
(4) 



The same definitions apply to the amplitudes for the CP 
conjugate cascade B + — > — > 7C + 7C~)D,ff + , with the 
change of weak phase sign — y in (4). We have set the 
strong phase of A B to zero by convention, so that 8 B is the 
difference of strong phases between the two amplitudes. 
For the CKM elements, the usual convention of the weak 
phases has been used. The value of \A B \ is known from the 
measurement of the B — ► D°K~ decay width using flavor 
specific decays of D°. The amplitude A(B~ — > D K~) is 
color suppressed and cannot be determined from experi- 
ment in this way [2]. The color suppression together with 
the experimental values of the ratio of the relevant CKM 
elements leads to the theoretical expectation r B <~ 0.1 — 0.2 
(see recent discussion in [7]). 

For the three-body D meson decay we define 

A D (s l2 ,s l3 )=A l2A3 e iS ^ 

= A(D Q ^K s ( Pl )n-(p2)n + (p 3 )) 

= A0 > ^K s (p 1 )7Z + ( P2 )7t-(p 3 )), 



(5) 



where Sjj = (pt + Pj) 2 , and pi,p 2 ,P3 are the momenta of 
the K s ,n~,7t + respectively. We also set the magnitude 
Ai2,i3 > 0, such that §12,13 can vary between and 2%. In 
the last equality the CP symmetry of the strong interaction 
together with the fact that the final state is a spin zero state 
has been used. With the above definitions, the amplitude 
for the cascade decay is 



A(B~ -» (K S X-X + ) D K-) = 

A B ^ D {A D (s n ,s n ) + r B e i{S *-^A D (j 13 , j 12 )) , 



(6) 



where is the D meson propagator. Next, we write down 
the expression for the reduced partial decay width 



dt(B- - {K S K-K + ) D K-) = (A 2 2il3 + r|A 2 3 )12 



+ 2r B Re A D (s l2 ,s l3 )A* D (s l3 ,si 2 )e 



-i(5 B -y) 



dp, 



(7) 



where dp denotes the phase space variables, and we used 
the extremely accurate narrow width approximation for the 
D meson propagator. 

2 In the following discussion we neglect D° —D° mixing, which is a good 
approximation in the context of the Standard Model, for more details see 
app. A of [8]. 




Figure 2. The partitions of Dalitz plot as discussed in text. The sym- 
metry axis is the dashed line. On the axes we have .S12 = mi 2 . % - and 

' S ' 13 = m K s n+ in GeV2 - 

The moduli of the D decay amplitude Ai 2i i3 can be mea- 
sured from the Dalitz plot of the D° — > K S K~K + decay. 
To perform this measurement the flavor of the decaying 
neutral D meson has to be tagged. This can be best 
achieved by using the charge of the soft pion in the decay 
D* + — ► D°7t + . However, the phase 5i2,i3 of the D meson 
decay amplitude is not measurable in B-factories without 
further model dependent assumptions. If the three-body 
decay D° — ► K s n~iz + is assumed to be resonance domi- 
nated, the Dalitz plot can be fit to a sum of Breit-Wigner 
functions, determining also the relative phases of the res- 
onant amplitudes. This is further discussed in section 4. 
The other option is to use data from charm factory, where 
weighted averages of the sine and the cosine of the relevant 
phase difference may be measured (see section 3). Here we 
assume that no charm factory data is available and develop 
the formalism without any model dependent assumptions. 

Obviously, to compare with the data, an integration over at 
least some part of the Dalitz plot has to be performed. We 
therefore partition the Dalitz plot into n bins and define 

ci = J [dp Ai2,i3Ai 3i i2 cos(<5i 2j i3 - <5i3,i 2 ), (8a) 

st = ^/>Ai2,i3Ai3 i i 2 sin(5i2 1 i3 - <5i3,i2), (8b) 

Ti = JdpA 2 l2l3 , (8c) 

where the integrals are done over the phase space of the 
i-th bin. The variables c, and contain differences of 
strong phases and are therefore unknowns in the analysis. 
The variables 7], on the other hand, can be measured from 
the flavor tagged D decays as discussed above, and are as- 
sumed to be known inputs into the analysis. 

Due to the symmetry of the interference term, it is con- 
venient to use pairs of bins that are placed symmetrically 
about the 12 <-» 13 line, as shown in Fig. 2. Consider an 
even, n — 2k, number of bins. The k bins lying below the 
symmetry axis are denoted by index i, while the remaining 
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bins are indexed with /. The /-th bin is obtained by mirror- 
ing the /-th bin over the axis of symmetry. The variables 
Ci,S{ of the /-th bin are related to the variables of the /-th 
bin by 



(9) 



while there is no relation between 7} and 7}. Note that had 
one used 12 <-> 13 symmetric bins centered on the symme- 
try axis, one would have had s, = 0. 

Together with the information available from the B + decay, 
we arrive at a set of Ak equations 



(10a) 



fr = Jdt{B~ -> {K s k-k + ) d K-) = 

Ti + r 2 B Tj + 2r B [cos(d B - y)a + sin(c5 B - y) Si ] , 

tj = j id t { B-^ { K s n-n- )D K-)= ^ 
Tj + r 2 B Ti + 2r B [cos(S B - y)c t - sin(<5 B - y)s,-] , 

f+ = Jdt{B+ -» {K s n-n + ) D K+) = 

T] + r 2 B Ti + 2r B [cos(5 B + y)c* - sin(5 B + y)s,-] , 

tf = [dt(B + - (K s %-% + ) D K + ) = 



(10c) 



(lOd) 



7] + r 2 B Tj + 2r B [cos(8 B + y)c/ + sin(<5 B + y)si] . 



These equations are related to each other through 12 <-» 
13 and/or y <-> — y exchanges. All in all, there are 2fc + 3 
unknowns in (10), 



q, s<, r B , 5b, y, 



(11) 



so that the 4k relations (10) are solvable for k > 2. In 
other words, a partition of the D meson Dalitz plot to four 
or more bins allows for the determination of y without 
hadronic uncertainties. This is our main result. 

When c, = or s,- = for all /, some equations become 
degenerate and y cannot be extracted. However, due to res- 
onances, we do not expect this to be the case. Degeneracy 
also occurs if 8 B = 0. In this case, y can still be extracted 
if some of the c; and/or s; are independently measured, as 
discussed in the following sections. 



3 Improved Measurement of c ( - and Si 

So far, we have used the B decay sample to obtain all the 
unknowns, including c, and s,, which are parameters of the 
charm system. We now show that the c, and s; can be in- 
dependently measured at a charm factory [10-12]. This is 
done by running the machine at the y/(3770) resonance, 
which decays into a DD pair. If one D meson is detected 
in a CP eigenstate decay mode, it tags the other D as an 
eigenstate of the opposite CP eigenvalue. The difference 




CP tag 



Figure 3. The interference between the decays of D mesons originating 
from ty/(3770) allow for a measurement of c; and s; at charm factories. 
Shown is a decay allowing for determination of c, . 

between the two decay widths gives [8] 
1 



Ci = - 



jdT{D% 



K s (pi)n (p 2 )7t + (p 3 )) 
K s (pi)n~(P2)n + (P3)) 



(12) 



where we have defined D|[ 



(£>° ± 25°) / V2. As stated 
above, obtaining these independent measurements reduces 
the error in the measurement of y by removing k of the 
2fc + 3 unknowns. 

In addition, if one of the D mesons decays into a non-CP 
eigenstate, we are sensitive to the s,- variables as well. Con- 
sider for instance a 1^(3770) decaying into afl,5 pair, of 
which one decays into KsTi + Ti~ and the other decays into 
some general state g. The partial decay width correspond- 
ing to the /— th bin of the K s n + n~ Dalitz plot and the j— th 
bin of the g final state's phase space is 

T lJ ^T l Tf + T I Tf-2(c l c^ + s^), (13) 

where T?, c S j, s g j are defined as in (8). In particular, if one 
chooses g = Ks% + ll~ and j = i (or j = i) one measures 
sf. If, on the other hand, g is a CP even (odd) eigenstate, 
s 8 j = 0, Tf = Tj = ±c 8 j and equation (13) reduces to (12). 

4 Assuming Breit-Wigner dependence 

If the functional dependence of both the moduli and the 
phases of the D° meson decay amplitudes Ad (sn, s i3) were 
known, then the analysis would be simplified. There would 
be only three variables, r B ,8 B , and y, that need to be fit 
to the reduced partial decay widths in Eq. (7). A plausi- 
ble assumption about their forms is that a significant part 
of the three-body D° — > KsK~K + decay proceeds via reso- 
nances. These include decay transitions of the form D° — > 
K s p° -> K s n-n + or D° -> K*~(S92)n + -» K s n-n+, as 
well as decays through higher resonances, e.g., /o(980), 
/ 2 (1270), / (1370) or ^(1430). An important feature 
is that there exists an overlap region of Cabibbo allowed 
D° -> K*-(892)7t + , D° -» K s p° decays, where the varia- 
tion of strong phase will be large, allowing for the extrac- 
tion of y. 

The decay amplitude can then be fit to a sum of Breit- 
Wigner functions and a constant term. Following the nota- 
tions of Ref. [13] we write 

A D (s 12 ,s 13 ) =A(D° ^K s ( Pl )n-{p2)n + (p 3 )) = 

= a e iS °+'Eare i5r ^ r (sn,su), (U) 
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where the first term corresponds to the non-resonant term 
and the second to the resonant contributions with r denot- 
ing a specific resonance. The functions stf r are products 
of Breit-Wigner functions and appropriate Legendre poly- 
nomials that account for the fact that D meson is a spin 
particle. Explicit expressions can be found in Ref. [13]. 

One of the strong phases 5, in the ansatz (14) can be put to 
zero, while others are fit to the experimental data together 
with the amplitudes a,. The obtained functional form of 
Au{s\2,s\^) can then be fed to Eq. (7), which is then fit to 
the Dalitz plot of the B ± — > {K s n^n + ) D K ± decay with r B , 
5b and / left as free parameters. 

5 Discussions 

The observables pr defined in (10) can be used to experi- 
mentally look for direct CP violation. Explicitly, 

a'^p = t~j - ft. = 4r B sin y [a sin 8 B T s ; cos 8 B ] , (15) 

Nonzero acp requires non-vanishing strong and weak 
phases. Due to the resonances, we expect the strong phase 
to be large. Therefore, it may be that direct CP violation 
can be established in this mode even before the full anal- 
ysis to measure y is conducted. With more data, 7 can be 
extracted assuming Breit-Wigner resonances (cf. section 
4). Eventually, a model independent extraction of 7 can be 
done (cf. section 2 and 3). 

The above proposed method for the model independent 
measurement of 7 involves a four-fold ambiguity in the ex- 
tracted value. The set of equations (10) is invariant under 
each of the two discrete transformations 

P a = {8B->8 B + n,y->y+n}, (16) 
P = {8 B - -8 B , 7^ -y,si - -st}. (17) 

The discrete transformation P n is a symmetry of the am- 
plitude (6) and is thus an irreducible uncertainty of the 
method. It can be lifted if the sign of either cos 8 B or sin 8 B 
is known. The ambiguity due to P can be resolved if the 
sign of sin 8 B is known or if the sign of 5; can be determined 
in at least some part of the Dalitz plot for instance by fitting 
a part of the Dalitz plot to Breit-Wigner functions. 

The rs suppression present in the scheme outlined above 
can be somewhat lifted if the cascade decay B — > DX~ — > 
(KsK~n + )DX~ is used [5,7]. Here X~ is a multibody 
hadronic state with an odd number of kaons (for instance 
K~n-%+,K-%° oxK s %-% Q ). The same formalism as out- 
lined above applies also to this case with trivial changes 
[8]. In addition to using different B modes, statistics 
may be increased by employing various D decay modes 
as well. An interesting possibility are Cabibbo allowed, 
D — ► K s n~n + n°, K~K + K S , and Cabibbo suppressed, 
D — > K~K + 7t°, n~n + n°, K s K + n~ decay modes, to 
which our formalism applies with very minor changes [8]. 



In conclusion, we have shown that the angle 7 can be de- 
termined from the cascade decays B — > {KgK~ k + )d- 
The reason for the applicability of the proposed method 
lies in the presence of resonances in the three-body D me- 
son decays that provide a necessary variation of both the 
phase and the magnitude of the decay amplitude across the 
phase space. The fact that no Cabibbo suppressed D decay 
amplitudes are used in the analysis is another advantage of 
the method and leads to a sensitivity on 7 at order 0(r B ). 
However, it does involve a Dalitz plot analysis with pos- 
sibly only parts of the Dalitz plot being practically useful 
for the extraction of 7. In reality, many methods have to 
be combined in order to achieve the required statistics for 
a precise determination of 7 [4]. 
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